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ABSTRACT 



Context. Magnetic fields are thought to be one of the possible mechanisms responsible for shaping the generally spherical outflow 
of evolved stars into often aspherical planetary nebulae. However, direct measurements of magnetic fields during the transition to the 
planetary nebula phase are rare. 

Aims. The aim of this project is to expand the number of magnetic field measurements of stars in the (proto-)planetary nebula phase 
and find if the magnetic field strength is sufficient to affect the stellar outflow. 

Methods. We used Very Long Baseline Array observations to measure the circular polarization due to the Zeeman splitting of 22 GHz 
H2O masers in the envelope of the proto-planetary nebula candidate star IRAS 19296+2227 and the planetary nebula K3-35. 
Results. A strong magnetic field of B\\ = -135 ± 28 is detected in the H2O maser region of the proto-planetary nebula candidate 
IRAS 19296+2227. The H 2 masers of K3-35 are too weak to detect circular polarization although we do present the measurements 
of weak linear polarization in those masers. 

Conclusions. The field measured in the masers of IRAS 19296+2227 is dynamically important and, if it is representative of the large 
scale field, is an important factor in driving the stellar mass loss and shaping the stellar outflow. 
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1. Introduction 

At the end of their evolution, a majority of stars go through 
a period of high mass loss that is an important source for re- 
plenishing interstellar space with processed materials. During 
the asymptotic giant branch (AGB) phase, this mass loss pro- 
duces circumstellar envelopes which are found to undergo a ma- 
jor modification during the rapid transition from AGB star to 
Planetary Nebula (PN). The standard assumption is that the ini- 
tial slow AGB mass loss in a short time changes into a fast 
superwin d generating shocks and acceleratin g the surrounding 
envelope (Kwok, Purton, & Fitzgerald, ll978l) . However, a large 
fraction of PNe have asymmetric shapes, with the majority of 
the young PNe being bipolar. Thus, at some point during the 
evolution to a PNe the AGB stars must undergo a process in 
which the spherically symmetric outflow is altered to produce a- 
spherical PNe morphologies. It has been shown that the energy 
contained in the outflows of young bipolar PNe is often orders 
of magnitude larger tha n can be provided by radiation pressure 
(Buiarrab al et all [200 lb . The source of this energy has been ar- 
gued to be magnetic fields, binary or disk interaction or a com- 
binati on of these (see e.g. Balick & Frank, l2002t iFrank et all 
2007, and refer ences therein). 

Maser polarization observations are the predominant source 
of information about the role of magnetic fields during the 
late stages of stellar evolution. Most observations have fo- 
cused on the masers in the envelopes of AGB stars, as OH, 
H2O and SiO masers are fairly common in these sources, and 
have revealed strong magnetic fields throughout the entire en- 



velope (e.g. lEtoka & Diamond, l2004t IVlemmings et all [2005; 
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Herpmeta]Jl2006). Because they are a short-lived (« 1000 yr) 
nansition phase between AGB star and PNe, proto-planetary 
nebulae (PPNe) are fairly rare. However, OH maser observa- 
tions reveal sim ilar strength magnetic field s as in their progenitor 
AGB stars (e.g iBains et al. L I2003L 120041) . Additionally, a very 
small fraction of the PPNe maser stars show highly collimated 
H2O maser jets. Currently 1 1 of these so-called 'water-fountain' 
sources have b een identifie d, and they are likely the progenitor 
of bipolar PNe dlmail 120071 and references therein). Polarization 
observations of the H2O masers of W43A hav e recently revealed 
that t he maser jet is magnetically collimated (IVlemmings et all 
120061) . 

Here we present the detection of a strong magnetic field 
in the H2O maser region of the circumstellar envelope of the 
PPNe candidate IRAS 19296+2227, only the second such de- 
tection in an object beyond the AGB phase. Additionally, we 
shortly describe the measure ment of linear polariz ation of the 
H 2 masers of the PNe K3-35 dMiranda etalll2001h , which was 
observed together with IRAS 19296+2227. The PNe K3-35 is 
only o ne of 4 PNe toward s which H2O masers have been de- 
tected (Suarezetall|2007). IRAS 19296+2227 is a source with 
both H2O and OH maser emission and IRAS colors indicative of 
an AGB star w ith strong m ass loss. Its near kinematic distance is 
estimated by Engelsj (I 19961) to be 3.2 + 1 .3 kpc, which implies an 
infrared luminosity of 5000 L.. Although this luminosity is con- 
sistent with that of an AGB star, its maser properties are atypical 
since both the OH and H2O spectra are narrow compared with 
the maser spectra of OH/IR stars. However, its maser properties 
do resemble those of a number of PPNe candidates. Therefore, 
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Fig. 1. Observed VLBA spectra as measured in the image plane of the H2O masers of IRAS 19296+2227 (left) and K3-35 (right). 



lEngelsl d 19961) concluded that IRAS 19296+2227 should also be 
considered a PPNe candidate. 



2. Observations 

The NRACQ Very Long Baseline Array (VLBA) was used on 
June 20 2003 to observe the H2O masers of the PNe K3-35 and 
PPNe candidate IRAS 19296+2227. The average beam width 
is * 1.0 x 0.3 mas and 1.0 x 0.7 mas at the frequency of the 
616 - 5 2 3 rotational transition of H 2 0, 22.235 GHz for IRAS 
19296+2227 and K3-35 respectively. We used 4 baseband filters 
of 1 MHz width, which were overlapped to get a velocity cov- 
erage of ta 44 km/s. However, for both sources maser emission 
was only detected in one of the observing bands. The data were 
correlated multiple times. The initial correlation was performed 
with modest (7.8 kHz= 0.1 km s" 1 ) spectral resolution, which 
enabled us to generate all 4 polarization combinations (RR, LL, 
RL and LR). Two additional correlator runs were performed with 
high spectral resolution (1.95 kHz= 0.027 km s _1 ) which there- 
fore only contained the two polarization combinations RR and 
LL. Such high resolution is needed to detect the circular polar- 
ization pattern of the H2O masers at sufficient accuracy. Each 
source-calibrator pair was observed for 6 hours. The calibrator 
was observed for 1.5 hours in a number of scans equally dis- 
tributed over the 6 hours. 



2.1. Calibration 

The data analysis path is described in detail in 
IVlernmings. Diamond. & van Langevelde! (|2002[) . It follows the 
method of iKemball. Dia mond. & Cotton (1995) and was per- 
formed in the Astronomical Image Processing Software package 
(AIPS). The calibration steps were performed on the data-set 
with modest spectral resolution. Delay and bandpass calibration 
were done using the calibrator sources J 1850+2825 and 3C84, 
while fringe fitting and self-calibration were performed on a 
strong maser feature. The calibration solutions were then copied 
and applied to the high spectral resolution data-set. Finally, cor- 
rections were made for instrumental feed polarization using the 
unpolarized calibrator 3C84. Unfortunately, no good calibrator 
data was available for the polarization angle calibration, making 
it impossible to properly determine the direction of the linear 
polarization vectors. 



1 The National Radio Astronomy Observatory (NRAO) is a facility of 
the National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 



The resulting noise in emission free channels of the modest 
spectral resolution I,Q, and U data cubes was ~ 4.0 mJy Beam" 1 
while that in the high spectral resolution I and V cubes was ~ 
8.2 mJy Beam 1 . In the channels with strong maser emission, 
the noise increased by ~ 10%. 

Fig.Q]shows the H 2 maser spectra for IRAS 19296+2227 
and K3-35 as measured in the image plane. The high-spectral 
resolution stokes I image cubes were analyzed using the AIPS 
task SAD, which was used to fit all maser features with peak flux 
densities higher than 8 times the rms in the channel map with 
two-dimensional Gaussian components. Subsequently, we only 
retained components that were detected in at least 10 consecutive 
channels, corresponding in total to ~ 0.27 km s _1 , as narrower 
features are unlikely to be real. To enable the most accurate po- 
sitiona l comparison with the observations in iMarvel & Boboltz 
(119991) we then also determined the mean right ascension and 
declination off-sets with respect to the reference feature using 
a flux density-squared weighting scheme. The peak flux, full- 
width half-maximum (Avi) and velocity ( Vlsr) were determined 
by fitting a Gaussian to the feature spectra. 



3. Results 

The result of our analysis of the H2O masers features around 
IRAS 19296+2227 and K3-35 is presented in Table. Q] Here 
we list the feature LSR velocity, full-width half-maximum (Av^) 
and the right ascension and declination positional offsets from 
the maser reference features (Aa and A6). A map of the maser 
distribution of IRAS 19296+2227 is shown in Fig. [3] Circular 
polarization of 1.1 + 0.2% was detected for the strongest maser 
(f) of IRAS 19296+2227 and its spectrum is shown in Fig.H No 
circular polarization was detected for the weaker features of K3- 
35. On the other hand, linear polarization at a level of 1 .6 + 0.2% 
was only detected on feature b of K3-35 while no significant 
linear polarization was detected for IRAS 19296+2227. 

The relation between circular polarization fraction Py, mag- 
netic field strength along the line of sight By and maser line 
width Avl is given by Py = 2Aff-B||/Avl- However, the coef- 
ficient App is a function of maser saturation level and for our 
analysis of the circular polarization of IRAS 19296+2227 / we 
us ed the full radiative transfer non-LTE inter pretation described 
in IvTemmings. Diamond . & van Langeveldel d2002l) to take this 
into account. In this interpretation, the coupled equations of 
state of the three dominant hyperfine compo nents of the 22 GHz 
H2O m aser transition was solved, following lNedoluha & Watsonl 
(119921) . for a linear maser in the presence of a magnetic field. 
It is then possible to fit the total intensity and circular polar- 



Table 1. Results 
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Fig. 2. Total power (/) and circular polarization ( V)-spectrum for 
the strongest maser feature (/) of IRAS 19296+2227. The thick 
solid line in the bottom panel shows the best model fit to V. This 
fit corresponds to a magnetic field strength of -135 mG along 
the maser line-of-sight. 

ization spectra directly to model spectra for different intrin- 
sic maser thermal line widths (Av t h). The emerging maser flux 
densities of the model spectra are expressed in 7/AQ, where 
Tj, is the brightness temperature (in K) and AO the beaming 
solid angle ( in sr). The uncer t ainties in these fits are further 
discussed in Vlemmings et al. (2006). The best fit model to 
the IRAS 19296+2227 spectrum gives By = -135 + 28 mG, 
Av th = 1.0 ± 0.3 km s'and log(7/ fc AQ) = 9.4 + 0.4. 

4. Discussion 

4.1. The H 2 maser shell of IRAS 19296+2227 

The distribution of the H2O maser features around 
IRAS 19296+2227 is shown in Fig. [3] The masers lie in 
an elongated structure with a length of ~ 40 mas. At a distance 



of 3.2 kpc, this corresponds to a linear size of ~ 130 AU. 
As no larger structure was found in VLA observations 
dMarvel & Boboltzl 119991) . this is likely the full extent of the 
H2O maser shell. 

iMarvel & Boboltzl d 19991) published the first VLBA images 
made at 3 epochs of observations in 1997 (May 31, June 17 
and July 10). Comparing our data with the earlier observations, 
we find that, while the flux has decreased by more than a fac- 
tor of two, the maser distribution is extremely similar. Our ob- 
servations fail to detect several of the weak features at Vlsr > 
41 km s~'but the two mo st prominent maser groups identified in 
IMarvel & Boboltzl ( 1 19991) have persisted for over 6 years. Still, 
the center velocity of some of the maser features, such as the 
strongest reference feature, has drifted by as 0.2 km s _1 , making 
it impossible to match up individual maser features. As a result, 
direct proper motion measurements cannot be m ade. If the maser 
motions detected by IMarvel & Boboltz] d 19991) are the genuine 
large scale motions, one expects, after 6 years of motion, a sig- 
nificantly different maser distribution. Fig. [3] indicates the posi- 
tions of the maser features for which proper motions were de- 
termined in IMarvel & Boboltzl ( fl9 99), transposed over 6 years. 
Assuming the maser reference feature was the same in all obser- 
vational epochs, the predicted maser distribution fails to match 
our observe d distribution, implying th at the proper motions de- 
termined in IMarvel & Boboltzl (1 19991) are likely dominated by 
random motions of the maser features. Assuming the maser fea- 
tures c and h in our observations originate from the maser clump 
located at -16.1 and 4.5 mas Rig ht Ascension and Declin ation 
off-set respectively in epoch 1 of IMarvel & Boboltzl dl999l) . the 
maximum proper motion of these features is ~ 0.15 mas yr~' 
which corresponds to ~ 2.3 km s -1 . 

As described above, the combined analysis of the total in- 
tensity and circular polarization spectra provides additional in- 
formation on the intrinsic conditions in the maser region. The 
best fit intrinsic thermal line width of the material that gives 
rise to the IRAS 19296+2227 FbO masers was found to be 
Av th = 1.0 km s" 1 . Since Av th « 0.5(7// 100) 1/2 , with T the 
temperature in the maser region, the H2O masers in the enve- 
lope of IRAS 19296+2227 exist at ~ 400 K. Additionally, we 
find that the emerging brightness temperature of the strongest 
maser feature is I^AQ « 2.5 x 10 9 K sr. Because the maser is 
marginally resolved we take the size of the maser to be 0.3 mas. 
The measured brightness temperature of IRAS 19296+2227/ 
is thus Tb ~ 2 x 10 11 K. This implies a beaming solid angle 
AD. « 1.2 x 10~ 2 sr, which is typical for the beaming angle o f 
the H2O masers around evolved stars dVlemmings et all 12005). 
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Fig. 3. The H2O maser region of IRAS 19296+2227. The large filled octagons are the maser features detected in our observations 
and are color-coded according to LSR velocity. The smaller circles with bl ack edges and vectors de note the predicted positions 
when transposing the maser features with proper motion determinations from Mar vel & Boboltzl (1 19991) between May 3 1 i997 and 
Jun 20 2003. The vectors indicate the maser feature trajector y, ending at the 1997 po sitions. This illustrates that the large scale 
maser proper motions are significantly less than determined in lMarvel & Boboita ([1999). 



Finally, as detailed in IVlemmings et alj (|2006), the emerging 
brightness temperature implies that the masers are unsaturated. 

4.2. Linear polarization 

No linear polarization was detected from the IRAS 19296+2227 
H2O masers down to ~ 0.4%. This is a similar limit 
as found for the Mira and supergiant stars in our previ- 



ous sample (V lemmin gs. Diamond. & van Langevelde! [2002; 



IVlemmings et al.L 12005). However, the strongest H2O maser of 
the PNe K3-35 was linearly polarized at a level of ~ 1.6%. 
Unfortunately, the lack of good polarization angle calibration 
made it impossible to determine the polarization vector orien- 
tation. The fractional polarization of the K3-35 maser is similar 
to that found for H2O masers in star forming regions and those 
in the 'water-fountain' PPNe source W43A IVlemmings et all 
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2006; IVlemmings & Diamond! [2006). Considering that the 
maser li near polarization fraction increases with maser satura- 
tion (e.g iDeguchi & WatsonL 11990 ). this seems to indicate that 
the saturation level of the H2O masers in the evolved stellar en- 
velopes is on average less than that of the maser in star forming 
regions, PNe and related 'water-fountain' sources. 

4.3. The magnetic field of IRAS 19296+2227 

The circular polarization of the strongest H2O maser fea- 
ture (/) of IRAS 19296+2227 implies a magnetic field of 
B|l = -135 + 28 mG in the H2O maser region. The 
field strength is similar to the fields previously detected us- 
ing H2O mase r polarization observations in the envelopes of 
evolved stars ( Vlemmings, Diam ond. & van Langeveldell2002t 
IVlemmings et aUl2005Ll2006l) . 

The origin of the magnetic field is unclear. As only one maser 
feature was strong enough to detect a magnetic field , a field lo- 
cal to the masers cannot be ruled out (e.g Soker, 2002). However, 
H2O , OH and SiO maser polarization observations of a number 
of other evolved stellar sources reveal ord e red large scale mag- 
netic fields (e.g. iBains et al L I2003L l2004t ICotton et al. L l2006t 
iKemball & Diamond 1 1 997t IVlemmings et all I2005L 12006). 

A field of 135 mG in a region with T « 400 K and a hydro- 
gen number density «h 2 of ~ 10 9 crrT 3 , typical for H2O masers, 
implies that the ratio between magnetic and thermal pressure 
B 2 KSnriYiJtT) =s 15. Assuming the masing material has a veloc- 
ity of w 3 km s , the ratio between magnetic and dynamic pres- 
sure B 2 /(8npv 1 ) ~ 5, where p and v are the density and velocity 
of the maser medium. Thus, if the measured field strength in the 
IRAS 19296+2227 H2O maser region reflects a larger scale field, 
the magnetic field is strong enough to play an important role in 
the dynamics of the stellar outflow and possibly also in driv- 
ing the mass-loss t hrough the effects of Alfven waves near the 
stellar surface (e.g.lFalceta-Goncalves & Jatenco-Pereiral 120021: 
IVidotto & Jatenco-Pereirall2006l) . 
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5. Conclusions 

We have measured a magnetic field of By = -135 + 28 mG 
in the H2O maser region of the PPNe candidate star IRAS 
19296+2227. The field strength is of the same order as the previ- 
ously measured field strengths in evolved stellar envelopes. This 
implies that there is no significant evolution of the magnetic field 
in the transition from AGB star to early PNe. The magnetic field, 
if representative of the large scale field, is strong enough to play 
an important role in the dynamics of the stellar outflow and pos- 
sibly also in driving the mass-loss. However a direct interpola- 
tion of the measured field strength to the stellar atmosphere is 
impossible. 

While no linearly polarized H2O maser emission is detected 
in the envelope of IRAS 19296+2227, ~ 1.6% fractional lin- 
ear polarization is found for the strongest maser feature of the 
PNe K3-35, implying that the masers aree at least partially satu- 
rated. With sufficiently sensitive observations, the measurement 
of linear polarization of the H2O masers of PNe will enable us to 
probe the magnetic field morphology early in the PNe phase. The 
circular polarization measurement and lack of linear polarization 
suggests the H2O masers of IRAS 19296+2227 are unsaturated. 
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